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Spherical,  Cylindrical  and  One -■dimensional  Plows 
of  Compressible  Fluids 

Joseph  3.  Keller 

I .  Introduction 

Relatively  few  boundary  value  problems  involving 
spherical  or  cylindrical  flows  of  compressible  fluids  have 
been  solved  exactly.   In  order  to  solve  more  problems  of  this 
tjrpe,  particularly  those  involving  variable  entropy,  we  have 
investigated  the  flow  differential  equations,  at  first 
without  regard  to  initial  or  boundary  conditions.   In  this 
way  a  class  of  non-isentropic  solutions  of  the  differential 
equations,  depending  upon  an  arbitrary  function,  has  been 
obtained.   Then,  in  the  second  phase  of  the  investigation, 
the  arbitrary  function  is  adjusted  to  satisfy  particular 
initial  or  boundary  conditions.   In  this  waj'"  the  free  expansion 
of  a  sphere  of  gas  into  a  vacuum  has  oeen   treated,  as  well  as 
the  -oropagation  of  finite  and  strong  shocks  in  variable  media. 
The  latter  treatment  includes  Frimakoff's  point-blast  solution 
as  a  sr^eclal  case. 

The  method  of  procedure  is  quite  simple.   In  section  II 
the  problem  is  formulated  in  Lagrangian  variables.   In 
section  III  a  class  of  solutions  is  obtained  by  the  method  of 
separation  of  variables.   In  section  IV  these  solutions  are 
specialized  to  isontropic  flow.   In  section  V  the  solutions 
are  applied  to  the  determination  of  strong  shock  waves 
propagating  in  variable  media.   In  section  VI  finite  shocks 
in  variable  media  are  considered. 

I I .  Formulation 

l-Je    consider  the  m.otion  of  an  inviscid,  non-heat- 
condiicting  fluid,  obeying  the  polytropic  equation  of  state. 
The  one,  two  and  three  dimensional  cases  will  be  treated 
together.   In  the  three  dimensional  case  y(h,t)  represents 
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the  radius  at  time  t  of  the  particle  vjith  the  Lagrangian 
coordinate  h,  which  is  defined  by  the  equation 

y(h,t) 

(1)  h  =  r       r'^"-p(r,t)  dr       n  =  1,2,3- 

Jy(o,t) 

In  the  above  equation  p(r,t)  is  the  density  at  time  t  and 
radius  r,  and  n  is  the  dimension,  which  is  3  ii"^  the  spherical 
case.   In  the  two  dimensional  case  y  represents  radial 
distance  from  an  axis,  and  in  one  dimension  y  is  a  cartesian 
coordinate . 

It  is  further  assumed  that  all  flow  variables  depend 
upon  y  and  t  onlj''  and  that  flow  occurs  only  in  the  y  direction. 
This  is  the  assumption  of  spherical,  cylindrical  or  planar 
syTETietry . 

From  the  definition  of  y,  the  particle  velocity  u  is 
given  by 

(2)  u  =  y,    . 

Similarly  from  (1)  the  density  p  or  specific  volujme  'Xl.   are 
given  by 

t  -W  -v_   -1    n-1 

^-5^  tT-  D   =  y   y^   . 

BecsiUse  of  our  assumptions  concerning  inviscidity  and 
non-conduction,  the  entropy  s  of  a  r^article  is  independent 
of  time  (at  least  between  successive  shocks).   Thus  we  have 

{h)  s  =  s(h) 

The  function  s{h)  is  given  by  initial  data  or  by  shock 
conditions,  and  is  assuiTied  to  be  knox^m. 

The  pressure  p  is  given  bj?-  the  equation  of  state 

(5)  p  =  P(p,s)  =  gCi^jS) 

For  a  polytropic  gas  or  liquid 


J . 
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(6)  g(t^,s)  -  g^  +  A(s)2^'Y   . 


The  xionction  A(s),  the  adiabatic  exponent  y»  ^^I'i  "^^'ic  Internal 
pressiire  e   are  assiomed  to  be  known. 

In  terms  of  the  above  defined  quantities,  the  equation 
of  motion  is 

(7)         y^t  =  -y'''^^&-J^7'"\)^  +  g^\] 

For  a  polytropic  gas  or  liquid  this  becomes,  using  (6), 

/Q\  ,/  \/  n-1   v-Y-l/  n-1   V    n-1   ,  ,  n-1   x-y   n-1 

(8)  y^^  =  YA(s){y   j^J  (y   y^^).^  y    '  ^-^ij       J^^'        Y 

Equation  (8)  is  a  second  order  partial  differential  equation 
for  5?-(h,t).   Tlie  coefficient  A  is  assumed  to  be  a  known 
f-onction  of  s,  and  by  [h^]  ,    of  h. 

The  problem  wc  consider  is  that  of  finding  solutions  of 
(8).   Prom  a  particular  solution,  the  flov;  variables  can  be 
foiand  hj   using  (2),  (3)»  iU-)    and  (5)' 

III.   Product  Solutions 

Let  us  seek  product  solutions  of  (3)  of  the  form 

(9)  y(b,t)  =  f{h)J{t) 

Inserting  (9)  into  (S),  and  separating  variables  wc  obtain 

(10)  j-    -Xj^^l-Y)-1  =  0 

(11)  -A[(f''-^f')-Y]'r^-2  -A'(f^-^f')"*^'  f""^    =  ^    . 

In  those  equations  X  is  an  arbitrary  separation  parameter  and 
primes  represent  derivatives  with  respect  to  t  in  (10)  and 
with  res'^ect  to  h  in  (11). 

To  solve  (10)  we  multiply  by  j'  and  integrate,  obtaining 

/To\  l/-'\2      \    .n(l-v)    1         /  ^  -,  \ 

(12)  ^(J  )   -  -^^.Yf   J  ^   ^'  =  ^a        (Y  ?^  1) 


\ 


V 


(13)  f(j')  -  X  lo-  j  =  ^a  (y  =  1)    . 


1....^   ,._..._  1 


Here  a  is  an  integration  constant.   Thus,  unless  j  is  constant, 
vfaich   is  possible  only  if  \  =  0  or  j  =  0,  we  have  the  solution 

/n)  N  r  r  2\    .n(i-v) ^  1'^^^,.    .    ,   /  n\ 

i^h)  J   ln(l-Y:!  ^  +H     aj  =  t    (y?^1) 

(15^)  J  [2\   log  j  +a]"-/^  dj  =  t         (y  =  1) 

The  integration  constant  in  (lli),  (Ip)  has  been  set  equal  to 
zero  vjithout  loss  of  generality  since  it  merely  shifts  the 
origin  of  time. 

To  solve  (11)  we  differentiate  out  and  obtain 

(16)  YA[f"f^"^4-(n-l)f^-^(f')^](f'^"^f')"Y"^  f'^' 

-  fn-2(fn-lf.)-Y  a'  .  ^    . 

We  now  consider  the  inverse  function  h  =  h(f)  and  denote 
h'(f)  by  q(f).   Then  (l6)  bscoraes 

(17)  YA[-q'q"-f'""^+  (  n-1)  r^"  V^  ]  [  f'^'^q"^  ]  "^"^  t""^ 

^  ^n-2r^n-l^-l]"  "^  A'(h)  =  \ 

We  shall  first  treat  the  case  y  t^   1  ^^^7   introducing  z(f)  and 
B(f)  defined  by 

(18)  z  =  q^"^   ,    B(f)  =  A[h(f)] 

•■•/ith  these  definitions,  (1?)  becomes  (prime  denoting 
differentiation  with  respect  to  f) 

(1-))   z'  H-  z[-(n-l)(Y-l)f"''"  +  ^^  (logB)'] 

+  ^(r-1)  ^(n-l)(Y-l)+l  ^0    ^ 
y3 


■■'  r 


■'-■.    ' .. 


'■■;■-, 


An  integrating  factor  of  (19)  is  f 
solution  is  (with  G  a  constant) 
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Y-l 
■(n-l)(Y-l)3  Y   ^^^  the 


(20) 


z  =  f 


(n-l)(Y-l)o  Y 


[G  - 


X(y- 


^r 


3 


-1/ 


Y 


d.f  ] 


Thus  we  have  q  from  (l8)  and  (20),  and  finally  since  h'  =  q, 


(21)   h  = 


J   fn-lg-l/Y|,._A(Y-l)  r  fg-l/Y  e^f 


-t-1 

r-1 


C 


df  +  E 


Here  E  is  an  integration  constant. 

Equation  (21)  gives  f(h)  implicitly,  and  thus  provides 
a  solution  of  (11)  for  y  ^  \ .      The  solution  mB.j   be  written 
more  simply  by  defining  ?(f5  by 


(22) 


P(f) 


V--1 


fB''^/Y  dfT'"' 
J 


if  \   f   0   this  can  be  solved  for  3(f)  and  yields 


(23) 


B(f )  =  (-\f}^(f  ) 


Y/f>'  >-Y  -.n 


The  solution  for  the  flow  variables  can  now  be  computed  from 

(2)-(5),  (9),  (21),  (22).  We  have  for  y  7^  1  and  \  ?^  0, 

remembering  that  f  =  yj~~  from  (9)j 

(21+) 


u(y,t)  =  yj'j"^ 


(25) 
(26) 


n-1 


.-1. 


ti(y,t)  =  -Xyr'"VF'(yj"^) 
?(y,t)  =  g^+  j"''^  F(yj"^) 


Equations  (2[i.)-(26)  give  the  flow  quantities.   In  fact,  these 
expressions  yield  a  solution  of  the  Eulerian  equations  of 
motion  for  an  arbitrary  function  F  provided  that  j(t)  is  given 
by  (lit-).   This  is  the  first  main  result  of  this  paper.   It 
is  to  be  noted  that  in  these  solutions,  u  is  proportional  to  y. 

In  the  excluded  cise  y  p^  1,  X  =  0  we  have  instead  of 
(2t^)-(26)  from  (2)-(5),  (9),  (21) 


n: 


J.,_' 
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(27) 
(28) 

(29) 


u(y,t) 
r(y,t) 


.T  .-1 

yj  .J 


-1 


=  yt 

1 


1  _' 


P(y,t)  =  g^  +  j-=^Y  G--Y 


=o 


g„*J!t-"v 


The  new  arbitrary  function  b  and  arbitrary  constant  S.    have 
been  Introduced  and  the  expressions  simplified  by  using  ( ll|-)  , 
which  gives 


(30) 


j(t)  =  a^'-t 


The  corresponding  solutions  when  y  -   \   are  obtained  from 
(17)  by  introdu.cing  3(f)  as  in  (l8),  after  vmlch  (17)  becomes 


(31) 

The  solution  of  (31)  is 


1-1       ,       ,  .  .-1    ,    B'         Xf         r. 

q    q         -(n-l)l         +   ™   +   —   =    0 


(32)  q(f)    =   r^"^3"-^(f)    exp  J       -Xf3"^(f)df      . 

Thus 


■X      r- 


(33)  h(f)    = 


f^-lB"^-(f) 


jxp   j       -\fB~^( 


f)df 


df   +  E 


Equation  (33)  yields  the  solution  implicitly.   Again 
we  define 


(3ij-)         F(f)  =  exp 
Then  if  X  7^  0  we  have 

(  35) 


/- 


XfB"-'-(f)df 


3(f)  =  -\fF(F')'- 


The  solution  for  the  flow  variables,  if  y  =  1  and  \   ^  0,    is 
again  given  by  (2[|.)-(26)  with  y   =  \   and  j(t)  given  by  (15). 
Similarly  for  y  =  1  find  X  =  0  the  solution  is  given  by 
(27) -(29)  with  y  =  1,  G  =  1  and  j(t)  given  by  ( 30) . 


1 


; : "  f: .. 


-r'       -.:,• 


!■ 


Equations  (2l4.)-(26)  and  (27}-(29)  represent  non- 
Isentropic  solutions  of  the  flow  equations  depending  upon  an 
arbitrary  function.   In  order  to  construct  these  solutions 
exp].icitly  one  need  merely  evaluate  the  integrals  in  ( Iti.) 
or  (15) •   Some  special  cases  of  these  integrals  yield  results 
which  are  listed  below  for  later  use. 

n(Y-l) 
2 

if  yj^l,    a=0,  11^  =  -1 
if  n(l-'^-)=  1,  ci,  p  arbitrary. 


if  y-n,    a=0,  iiA-p-i  /  -1 


(37)   3{t)=e' 


+aV2, 


:X 


(38)   j(t)  =  |t^+  at  +  B 


IV .   Isentropic  Case 

The  above  solution  can  be  specialized  to  the  isentropic 
case  by  setting  A  =  B  =  constant.  Then  frora  (22)  and  ( 3L|.)  we 
obtain 

(39)   F(f)  -   G  - 


A(v-1).2]Y-1 


(IiO)   F(f)  =  exp 


2y3 

-\f^ 
23 


V7 


Y  7^   1 


Y  =  1 


Thus  the  solutions  given  "oj    (21!.) -(26)  become,  for  y  7^  1  ^^^ 
\   =   0 


(ill)  u(y,t)  =  yj'j"^ 


([^2)  r(y,t)  =  jV/^Tr.  .  MjzJ-)   y2j-2 


2y3 


1/y 


-1 


Y-1 


(,Il3)  P(y,t)  =  g^+  j 


-ny 


G  -  ^W-1)  ,2  -2 
2y3^/y  " 


V 

V-1 


In  these  solutions  B  and  G  are  constants  and  j  is  given  by  (lii). 
For  \  =   0,    and  all  y  the  solution   (27)-(29)  applies  with  b 
constant. 


■■Jl 


•  I.   -    ;   I 


■'.    •.  ! 
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For  Y  =  1  ^^^  ^^^  ^j  ('-i-1)  is  unchanged,  but  the  other 
equations  become 

^  .-2  2 
(Ml.)   r'(y,t)  =  j^B  exp  -^^.^/ 


-2  2 


ihS)      P{y,t)  =  g^+  j-^  exp  2Li^ 


23 


Here  B  is  a  constant  and  j  is  given  by  (15). 

As  an  example  of  these  isentropic  solutions,  let  j(t) 
be  given  by  (36)  and  assuraei  G  =  0  in  (1^-2),  (14-3).   Then  (l+l)' 
(I.I.3)  yield  the  power  solutions  (if  y  7^  1,  --^JCI — L  ^   _i) 


(i;6)   u(y,t)  =  ^TITFYTT?  yt"^ 

1         -2 

(ii7)  t;(y,t)  =  [i^(ll!.-ir^l)2]^r-i  (yt-i)^'-i 


(i;8)   p(y,t)  =  g^+  B 


rB,n(Y-l;+2.2|l-Y  /^.t-^)^'"^ 


V.   Application;   Strong  Shocks  In  Variable  Media 

Let  ixs  suppose  that  a  shock  given  by  the  equation 
y  =  ?. (t)  moves  into  a  variable  medium  of  density  p  (y), 
pressure  P^^ly)  3-nd  velocity  zero.  '2''ae   pressure  p,  velocity  u 
and  density  p  just  behind  the  shock  are  related  to  the 
corresponding  quantities  in  front  of  it  by  the  shock 
conditions.   These  conditions  are 


(y+i)p+  (^-i)p, 


(i+9)     f~  =  TFTTp-+T7fTT 


p. 


Y+1 


(y  ^  -1) 


(50)      u  = 


2(P   -P^) 


r        "5 


£2p^[(Y+l)p+  (Y-l)Pol} 


~TT7o 


p 


-11/2 


Ty+i.'p 


oj 


(  51)         R 


_  r(Y^i)PMY-i)P,ii/2^  r(-,.i)r.ii/' 


L 


2P, 


O": 


.;  ( .  L  ."V 
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The  second  expressionfB  on  the  right  apply  to  a  strong  shock, 
for  T>;hich  p  »  p^. 

We  shall  now  assLijr.o  that  the  flow  behind  the  shock  is 
a  "product"  solution  given  by  {2i|)-(26),  and  that  the  shock 
is  strong.   V/e  wish  to  determine  the  functions  P(yj~  ),  R(t), 
D  (y)  and  p  (y)  for  which  such  a  solution  is  possible.   To 

'00 

this  end,  we  insert  (21)-) -(26)  into  (ij.9)-(5l)  and  obtain 


/  i-'o  ■' 

w2; 


(53) 


^2g^-.2J-^Yp(Rj-l)Jl/2 


"(^+iiTjrRT 


i5k) 


R  - 


(Y+DiSo^  r^^F(Rj-i)}' 


2p^(R]- 


1/2 


Prom  {S3)    ^^nd  (5^-!-)  wc  have 


(55) 


R   "■  2    j 


Thus 
(56) 


R(t)  =  RQ[j(t)] 


^ 


Inserting  (56)  into  (5.3)  and  making  use  of  (12), 
satisfied  bj   j(t),  we  find 


(57)   F(x) 


^P. 


fr-i   ^i-Y^ 


y 


•[HrirTy^-^Cr^'-Dir 


-2nY   (2n+2)Y-2 


v-1 

R'     X 


2n 


-«o(-%'''^"'- 


Thus  R  is  determined,  P  is  related  to  p  ,  and  (53) >  (5^)  are 
satisfied,  although  R  ,  D,  and  X  are  still  arbitrary.   We  now 
Insert  (56),  (57)  into  (52)  to  determine  P.   Me  obtain,  if 
go  =  0. 


r  r--^  ■•■    I  ;        r.     O  ' 


o\:      .  ■  -L . 


-L    \ 


I"  ^ 
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(58)      P(x)    =  F^x'^ 


•    2XR2n  -1-1/2 

ni l-y) 


Prom    (^y)    ^nd    (58)    we    then   find 

2nv 


(59)      D^(R) 


2F_R 


Y-1    - 


o"  o 


2\R 


2n 


Y+1        [n(l-Y) 


+  brT'*'"'"  R 


l\n     2n(Y-l)  n  -3/2 
Y+1 


( n-2)Y-3nH-2      2(  n-2 )  y-Gn+l,^ 


.R 


Y+1 


R. 


(  Y+1 )( Y-1) 


In  (58)  and  (59)  F   is  an  arbitrary  constant.   The  corres- 
ponding results  with  g   ^  0  are  somevrhat  more  comDlicated. 
We  have  thus  obtained  a  solution  with  a  strong  shock 
moving  into  a  variable  medium  at  rest,  v;ith  density  given  by 
(59)-   The  constants  F    ,    R  .  \  and  D  are  arbitrary  in  this 
equation,  but  only  tv/o  essential  combinations  of  these 
constants  occur.   The  shock  curve  is  given  by  (56),  and  j(t) 
by  (III.).   The  flow  is  given  by  (2)4.) -(26)  with  F(x)  given  by 
(58).   The  flow  might  be  produced  by  a  piston  following  one 
of  the  particle  raths 


(60; 


Y(t)  =  Y^j(t) 


The  solution  was  deduced  for  y  /  -1  ^^'^  ^  /  ^'« 

As  an  example  of  these  solutions,  lot  us  suppose  that 
p  ( R)  =  constant.   From  (59)  we  find  that  this  requires 


(61) 


D  =  0 


n-2 


The  last  condition  can  be  fulfilled  only  for  n  =  3,  in  vjhich 
case  Y  =  7.   (n  =  1  leads  to  v  =  -1,  v:hich  V7as  excluded  in 
the  derivation,  and  n  =  2  yields  no  value  of  Y' )   In  this 
case  we  have  from  (1)-|-),  (58) 


(62) 


J  = 


L  3   ^J 


Tl/10 


-1/  \    3   00 
F(x)  =  x-^  — — 


•r  -    r     •.•   I  - 


1-    •  :■  •  ■ 


11. 


The  solution  computed,  by  using  (62)  in  (21j.)-(26)  (lAfith.  X  <  0) 
is  exactly  the  point  blast  solution  of  Taylor's  type  first 
found  by  H .  Primakoff,  and  applicable  to  high  energy 
explosions  in  water.   [see  [I],  p.  h^h] 

VI .   Application;   Finite  Shocks  In  Variable  Media 

He  may  apply  the  above  method  to  determine  finite 
shocks  in  variable  media.   Then  we  must  satisfy  the  exact 
shock  conditions  (Ij.9)-(51)  rather  than  the  strong  shock 
conditions.   Inserting  (2k)-(26)  into  {Ll9)-(51)  yields  with 

g-       ~    O, 

(63)   — L-^_  = °   ,   X  =  Pj  ^- 


i^k)  Rj  J 


^"'"^-Po 


?[(Y+1)J"'^F+(Y-1)PJ^^^ 


(65)         R  =  -1-  [(v+l)j-^Yp+ (y-Dp^]^/^   . 

J2Po 

Equations  (63) -(65)  are  a  set  of  two  first  order 
ordinary  differential  equations  and  one  algebraic  equation 
for  the  determination  of  the  four  functions  F(x) ,  R(t),  p  (R) 
and  p  (R).   Equation  (I/4-)  gives  j(t).   This  system  is  evidently 
ujider determined,  and  will  have  infinitely  many  solutions.   In 
the  strong  shock,  case,  however,  p   did  not  occvir  and  thus  the 
system  was  determined. 

To  solve  the  above  system  we  could  impose  some  relation 
between  p  and  p   and  eliminate  both  these  functions  by 
mea.ns  of  that  relation  and  {Sl\.)  .      A   pair  of  simultaneous 
first  order  equations  for  P(x)  and  R(t)  would  result.   However 
these  equations  can  be  treated  separately  since,  from  (6L|.) 
and  (65)  wc  have  by  eliminating  F, 


'J  ^ 


\       :     ^       ~     ' 
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(66) 


^  =  ^ 


T5  -: 


I  .-1 


H 


b:±^(Hj'r^-)^^Y.oo;^, 


1/2 


-1 


Th5-S  is  an  equation  for  R(t)  if  the  ratio  PqO    is  given  as 

a  function  of  R.   (Equation  (66)  also  holds  when  g  /  0.) 

After  sol\'-ing  this,  (63)  can  be  solved  for  P,  and  then  p 

and  p   can  be  obtained. 
'^o 

We  v/ill  now  investirate  the  shock  curve  R(t)  given  by 

(66),  in  the  special  case  in  which  j(t)  is  given  by  (36)  with 

_n     2 
C  =  0  ana  "rPoP  "  ~  "^o  ^^  constant.   The  latter  assumption 

means  the  temperature  ahead  of  the  shock  is  constant.   Making 
use  of  (36),  (66)  becomes 


(^7)   H  =^rH^^,— Rt-1  ± 


>2 


IXLil 


U[n(Y-l)+2] 


-1,2     2 
■2   (^^^  ^''   ^  -o 


1/2 


If  we  nov;  introduce  '''T(t)  -  Rt'   in  (67)  we  obtain  an 
equation  for  U  in  which  the  variables  separate.   One  solution 
is  the  constant  solution 


(63) 


U  =  1  c  (l--.l±^)-l/2 
o^    n( y-l )+2 ' 


In  this  case,  which  is  physically  possible  only  when 

%•  + 1 
—7 — '-ti — ?r  <  1,  tne  shock  curve  is  the  straip;ht  line 
n(v-l)+2 


(69) 


R  =  IJt  =  ±  c^(1--7-14-Tt)'"^^^  t 
o'    n( y-l ; +2' 


Vfhen  n   does   not   have    the    value    e-iven  by    (68),    the 
equation   for  U   yields 


(70) 


/(t) 
j  |_^2H(Y-l)+i4 


1)U 


(v+l)' 


[2n(Y-l)+l|] 


P   IT^  +  c^ 
2  oj 


1/ 


dU   =   log 


Integrating   (70)    yields 


13. 

1-a   .       -1   ,     L 

(71)   bt  =  I  1  a  +  (a-l)cos 


Here  0  =  cot   -^  ,  a  =  p^>  ^_y  ^  j  and  b  is  a  constant.   Thus  we 
have 

(72)  H  =  TTt  =  -^  .  bt  cot  © 

Since  bt  is  e;iven  in  terms  of  9  by  (71)>  H(  t)  is  given 

ab 
parametricaily  in  terms  of  0.   A  p:raph  of  — R  versus  bt  is 

given  in  Fip'ure  1  for  n  =  3  ^nd  y  -    l.I}.,  which  represents  a 
spherical  shock  in  air.   When  the  minus  sign  is  chosen  in 
(^7) -(71)  the  same  curves  are  obtained  as  with  the  plus  sign, 
provided  t  is  replaced  bj/  minus  t. 
For  small  values  of  t  we  have 

(73)  R  »  Kt'^^ 

where  K  is  a  constant.   For  y  =  7  '^^nd  n  =  ?  we  find  a  =  -^   and 

thus  R  ^  X   ^     ,    which  is  the  behavior  of  Frimakof f ' s  point 

blast  solution  for  small  t.  For  large  t,  R  behaves  linearly 

in  t ,  as  in  ( 69 ) • 

To  determine  p   and  F,  we  have  from  (65)  and  the  condition 
_i     2 
YPopQ^  =  ^o' 

•  2 
Sd^R  -{y-Dp^     2yP   .n  1 

(?U)      F  =  ^  \^^^-   =   -_I^  R2jnr  _  iz^    .nr    . 

(^+i),]    ^  VY+1)C-  ^^^ 

Thus,  if  dot  denotes  time  derivative  and  prime  denotes 
derivative  of  a  function  with  respect  to  its  argument,  we  have 

(75)   P  =  _^:^(p;H3.nY^2p^^^.nY^^^p^^2,nY-l..j 
(y+1)Cq 


Y- 


Y+ 


i(p;Rf^-nYP,j^^-^j') 


'  •  -1     -P  t 
=  F  (Hj  -"  -Rj  "^j  ) 


TTsing   (7U.),    (75)    in    (63)    yield?    the    following   equation   for  p^ 


, '  ■    ( 


111. 


(76) 


Y+l  "  o 


.ny-l ..  t  \  I 


o 


L(y+1)c 
^(P^R.I  ^  +  nYPoJ  '   J  )j  • 


-XRj'^^p^c-^lRj-^-  Rj^i') 


T-1 


o  o 


°    L       (Y-M)cr    ^'"' 


-1 


Upon  introducing  p   =  p'R,  the  above  equation  can  be  integrated 
yielding  p  [R{t)]  as  a  function  of  t .   To  obtain  p  (R),  we 
must  insert  t  =  t(R)  from  (71),  (72)  into  this  result.   Then 
i-jhen  p  [R(t)]  is  known,  F  can  be  found  from  (7I4-).   These 
calculations  v;ill  not  be  carried  out  here. 


VTI.   Free  Expansion  and  Other  Problems 

Tn  this  section  some  flows  described  by  the  preceding 
"product"  solutions  will  be  discussed.   First,  let  us  consider 
the  flow  given  by  (Ij.1 ) -( L(.3)  with  g   =0.   This  describes  the 
isentropic  flow  of  a  gas.   If  .]*'(0)  =  0,  it  starts  from  rest. 
If  the  constant  G  is  chosen  so  that  p  =  C  for  some  value  of 
yj"  ,  then  p  remains  zero  for  this  value  of  yj   ,  which 
therefore  represents  a  "free"  surface  of  the  gas.   Thus  the 
flow  represents  the  free  expansion  of  an  isentropic  gas  into 

y 

vacuum.   The  equation  of  the  free  surface  is  y  =  . ,  9 'v  j(t), 

where  v   represent ?:  the  location  of  the  free  surface  at  t  =  0. 
"  o 

In  a  similar  way  (2[|.)-(26)  may  represent  the  free  expansion 
into  vacuum  of  a  gas  xvith  initially  variable  entropy,  provided 
g   =  0  and  F(yj~  )  vanishes  for  one  value  of  its  argument. 
For  certain  choices  of  j(t),  this  solution  can  also  represent 
the  collapse  or  expansion  of  a  hole,  or  the  initial 
contraction  of  a  sphere  of  gas  toward  the  origin  due  to  an 
initial  inward  motion,  followed  by  momentary  rest  and 
subsequent  expansion.   When  g   is  not  zero,  similar  motions 
of  a  liouid  are  described. 


r .  r- 


■j  ■'■ 


cv 


i." 


I  t  : 
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VTTI.   Other  Solutions 

In  the  preceding  sections  procuct  solutions  of  (3) 
were  considered.   In  this  section  tv;o  other  kinds  of  solutions 
will  be  examined.   First,  let  us  assume  that 

(77)  y  =  f( at  +  Sh) 

Here  o  and  f?  are  constants.   We  find  for  the  flow  quantities 


(73)    ^.       u  -  y,  =  cf  =  cf[f  "(y)] 


(79)  t;  =  p   =  Be"  y  ~  u 

(30)      P  =  Ap^^  =  A(aP"^y^"%"^)Y 

We  see  that  the  flow  quantities  are  independent  of  t ,  or  in 
other  words,  that  a  solution  of  the  form  (77)  yields  steady 
flows . 

Inserting  (77)  into  (S)  yields 

(81)   a^f  =  YA(3f^-^f)-^'-^(6f^-^f)>.f^"-^  +A^(Sf'"-^n-Yf'^-^   . 

Since  f  is  a  function  of  at  + Bh  and  t  doesn't  appear  in  (3l), 

h  cannot  aij-Dcar  either.   Thus  A  is  constant  and  A,  =  0.   Now 

"  ■  h 

if  X  denotes  the  argument  of  f,  and  if  we  consider  the  inverse 
function  x(f),  vie    obtain  from  (Si)  the  following  equation 
for  q(f)  =  x(f) 

(32)   -c^qq--^  =  YA(Pf^'-^q'^)"'^-^(-6f^-^qq-'+B(n-l)f^-2q-2)pfn-l. 
This  eauation  is  exact  ri.nd   yields  upon  integration 

(83)    bf(l-v)(-^)aJ--4q-2  .D  ,   b=XAP^  ,   ^^1 

0  (y-l) 

(3J|)    %-  q  ^  +  A  log  q  =  (n-l)Af  +D         .       ,    y   =   1       . 
Here  D  is  a  constant. 


;■■  .■  ir. . 


I       ,  -  t, 


'    .-.  :^■      :;■ 


.1  .■•■:J 


:   "I  . 
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From  (78)  r.nc  the  facts  that  f  =  q"  ,  f  =  y,  we  have 


(85) 


u 


qTyT   • 


Thus  the  steady  solutions  are  obtained  from  the  algebraic 
equation  (83)  or  (3[|.). 

A  more  Interesting  type  of  solution  is  given  by 


(86) 


y  =  t^f(K(h)t) 


Here  a  is  a  constant,  and  f,K  are  functions  to  be  determined, 
These  solutions  are  like  the  tyoe  considered  by  Guderly, 
Taylor,  von  Neumann  and  Calkin  and  Courant  and  Priedrichs. 
Inserting  (26)  into  (3)  yields 

(37)   (a-l)at^"^f  +  2aKt"'-'-f  +  K^t"*f 


YAlt^^^lf^-lKf 


-Y-1 


, a(n-l) pn-l, an+1 


f   ^fK+  f  ^YTti  +   (n-l)K  f  tf" 


,a(n-l) 


in+1  ^n-1-',  ?,  I  ~Y 


+  A,,t--^--f'^-l[t-^^^-^f'^-^KfJ 


In  order  that  h  occur  in  this  equation  only  in  the  combination 
X  =  p[(h)t,  we  find  the  follo;\ring  expressions  for  K(h)  and  A(h) 


(35) 


K  =  \pO-+  C  j^'^ 


A  =  DK'  "K 


Y-Vo  +  1 


Here  a  =  an(  l-y) -y-'-'i+S ,  D  is  a  constant  and  A  is  constant  if 
e( Y-1) +0+1  =  C. 

Under  these  conditions,  (67)  becomes 

(39)   (a-l)af  +  2axf  +  x2f  =  xS^^^"^^  ^  ^-"^^f-Vc 

r  Z'.  ^  '^'^  ^   (n-l)xf  A    ,   ,  ,      J 
1^-^  (  e  +  -^  +  -^ — -^ 1  +  e  (  Y- 1 )  +  c  +  1  j 


■M.  . 


■v:.     ( 


':    If 


17^ 


To   reduce    (39)    to   a   first    order   equation,    we    follow 

von  Neumann   and    Calkin,    by   letting   5    =    (  a+^y )  [  2  +  n(  y-l )  ]         and 


(90) 


X   =   e 


e      F 


Tcie  n   (39)    be  c  ome  s 


(91)       (a-l)aP+  2a(-^  +5)F+  (;^  +  5-l)(4-  +  5)P 


as 


ds 


QS 


r  .  1-Y 


=  P(n-l)(l-r)       d    ^ 

LOS 


ye  + 


v(n-l)(|^+5)F 
+ ^ +   eY~e  +  c+l 


If  we  now  let  q  =  (-r-+o)F  we  oDtain  a  first  order 
equation  for  q  in  terms  of  F. 

(92)    [l-  YrF(^-'^^^l-^^q-'^-^J(q-cF)4 

=  P'^"-^' '  ^'^■'q~'^D[Yf:  +  Y(--l)qF"-^+eY+a-i-e]-(a-l)aF-2aq  . 


Special  cases  of  this  equation  have  Deen  studied  by 
von  Neumann  and  Calkin. 
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